Body-wave interferometry using multi-dimensional deconvolution after wavefield decomposition at free surface

Passive seismic methods using earthquakes can be applied for extracting body waves and obtaining information of subsurface structure. To analyze earthquake data, one traditionally uses receiver functions, tomography, and seismic interferometry. In this study, we retrieve reflected body waves by applying seismic interferometry to the recorded ground motion from a cluster of earthquakes. We employ upgoing/downgoing P/S wavefield decomposition, time windows, time reversal, and multi-dimensional deconvolution to improve the quality of the extraction of reflected body waves with seismic interferometry. The wavefield separation and seismic interferometry based on multi-dimensional deconvolution (MDD) allow us to reconstruct PP, PS, SP, and SS reflected waves without unwanted crosstalk between P and S waves. From earthquake data, we obtain PP, PS, and SS reflected plane waves that reflect off the same reflector, and estimate Pand S-wave velocities.